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We have studied the sound perturbation of Unruh's acoustic geometry and we present an ex- 
act expression for the quasinormal modes of this geometry. We are obtain that the quasinormal 
frequencies are pure-imaginary, that give a purely damped modes. 
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I. INTRODUCTION 

m : 

General relativity has many successful predictions, some of them are passing experimental tests with high precision, 

■ for example mercury's perihelion and gravitational redshift. On the other hand, the most spectacular and promising 
prediction is the existence of black holes. In this sense the theoretical area of black hole physics became an attractive 
and productive field of research in the past century. Those classical and quantum properties are well understood 
within the general relativity framework. From a classical point of view, nothing can escape from the black holes; 
however, from the quantum point of view, black holes are not completely black and they can emit radiation with a 
temperature given by h/SirksGM 0. This radiation is essentially thermal and they slowly evaporate by emitting 

\ quanta. Another property of interest that play an important role in black hole physics are the quasinormal modes 

■ (QNM's). They determine the late-time evolution of fields in the exterior of the black hole, and numerical simulations 
of stellar collapse and black hole collisions have shown that in the final stage of such processes the quasinormal modes 
eventually dominate the black hole response to any kind of perturbation. QNM's has received great attention in 
the last few years, since it is believed that these models shed light on the solution or understanding to fundamental 
problems in loop quantum gravity. In particular, from Hod's proposal 0, modes with high damping have received 

jy-^ great attention, specially their relevance in the quantization of the area of black holes 0, 0, 13 and the possibilities of 
I fixing the Immirzi parameter in loop quantum gravity 0, 0, Q ■ 

Obviously the experimental test are impossible to make at the level of terrestrial laboratories. This opens the 
interest in analog models that mimic the properties of black hole physics. In this sense, the field of analog models of 
gravity allows in principle that the most processes of black hole physics can be studied in a laboratory. In particular 
^jjy the use of supersonic acoustics flows as an analogy to gravitating systems was for first time proposed by Unruh [t| 
and with the works of Visser 0, EH E3 EJ El has received an exponentially growing attention. 

The basis of the analogy between gravitational black holes and sonic black holes comes from considering the 
propagation of acoustic disturbances on a barotropic, inviscid, inhomogeneous and irrotational (at least locally) fluid 
flow. It is well known that the equation of motion for this acoustic disturbance (described by its velocity potential ) is 
identical to the Klein-Gordon equation for a massless scalar field minimally coupled to gravity in a curved spacetime 
|l2l liH IT3|. QNM's modes for acoustic black holes in 2+1 dimension was computed in Refs. [TEl flU IT7L ITq and for 
analogous black holes in Bose-Einstein condensate in Ref. [l9|. In this work we analytically compute the QNMs or 
acoustic disturbances in Unruh's 3+1 dimensional sonic black hole (Laval nozzle fluid flow), and also we show their 
large damping limit. 

The organization of the paper is as follows: In Sec. II we specify the Unruh's Sonic Black Hole. In Sec. Ill we 
determine the QNMs, and the large damping limits. Finally, we conclude in Sec. V. 



II. UNRUH'S SONIC BLACK HOLE 



The idea of using supersonic acoustic flows as analog systems to mimic some properties of black hole physics was 
proposed for the first time by Unruh |9j- Essentially he showed the possibility to use sonic flow in order to explore 
properties like the Hawking temperature near the sonic horizon and, in principle to developed an experimental 
setting to study the fundamental problem of evaporation of real general relativity black holes. As we mentioned in the 
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introduction, the basis of the analogy between gravitational black hole and sonic black holes comes from considering 
the propagation of acoustic disturbances on a barotropic, inviscid, inhomogcncous and irrotational (at least locally) 
fluid flow. It is well known that the equation of motion for this acoustic disturbance (described by its velocity 
potential ■0) is identical to the Klein-Gordon equation for a massless scalar field minimally coupled to gravity in a 
curved space 0, 0, Q , 

In Unruh's work, the acoustic geometry is described by the following sonic line element 



ds 2 = f (- (? - «?) dr 2 + + rW) , (1) 

where po is the density of the fluid, c is the velocity of sound in the fluid (by simplicity we will assume these quantities 
constant) and Vq represent the radial component of the flow velocity. 

On the other hand, if then we assume that at some value of r = r+ we have the background fluid smoothly exceeding 
the velocity of sound, 

Vo = -c + a(r-r + ) + tf(r-r + ) 2 , (2) 

the above metric assumes just the form it has for a Schwarzschild metric near the horizon. In this limit metric our 
reads as follows 



ds 2 = El ( _2aS(r - r + )dr 2 + f 7 * r + r 2 dil 2 ) , (3) 

c \ 2ac(r — r + ) J 

where a is a parameter associated with the velocity of the fluid defined as (V • ~v)\ r =r + H3- Note that this geometry 
was study in Ref. 20j where the author studied the low energy dynamics and obtained the greybody factors for the 
sonic horizon from the absorption and the reflection coefficients. 

In the quantum version of this system we can hope that the acoustic black hole emits " acoustic Hawking radiation" . 
This effect coming from the horizon of events is a pure kinematical effect that occurs in any Lorenzian geometry 
independent of its dynamical content |l3j . It is well known that the acoustic metric does not satisfied the Einstein 
equations, due to the fact that the background fluid motion is governed by the continuity and the Euler equations. 
As a consequence of this fact, one should expect that the thermodynamic description of the acoustic black hole is 
ill defined. However, this powerful analogy between black hole physics and acoustic geometry admit to extend the 
study of many physical quantities associated to black holes, such as quasinormal modes which we consider in the next 
section. 



III. QUASINORMAL MODES OF UNRUH'S SONIC BLACK HOLE 

The basis of the analogy between Einstein black holes and sonic black holes comes from considering the propagation 
of acoustic disturbances on a barotropic, inviscid, inhomogeneous and irrotational fluid flow. It is well known that 
the equation of motion for these acoustic disturbances (described by its velocity potential <I>) is identical to the Klein- 
Gordon (KG) equation for a massless scalar field minimally coupled to gravity in a curved spacetime 0, 0, 0] , 



(vW^$) = 0. (4) 

In order to compute the QNM's we apply the standard procedure described in Refs.^30|[23' an d begin to rewrite 
the metric Q in the following form 



ds 2 = -f(r)dr 2 + ^-+r 2 dn 2 1 (5) 

where f(r) = 2ac(r — r + ) In order to compute the QNM'S we need to solve the Klein Gordon equation Q in curve 
space described by the metric (|SJ) and, by virtue of symmetries of the metric we use the following Ansatz for the scalar 
field 
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$ = e- iu)t Y l m {6,ip)R{r), 
with this ansatz the KG equation can be separated as follows 



-J- (d e {smedeYne^)) + -^tfxrV,?) = Kl + W m {0,v), (6) 
sinf \ sine/ * J 

1 d (r 2 f{r) d R{r)] + ( bj2 + ^ + ^\ R(r) - 
c 2 r 2 dr \ dr J \f(r) r 2 J 

Now, if then we consider the change of variables z = 1 — r + /r, that transform the radial equation to 

, d 2 R dR „, . _ , . 

z(l-z)— + — +P(z)R = 0, (7) 
dz z dz 

here, 

P(z) = ^—,-A (8) 

2(1 - Z) 

where, 

A = (9) 
Zar+ 



B = — 

1 2a 

Note that in the new coordinate system, z = correspond to the horizon and z = 1 corresponds to infinity, and if 
we consider the definition 

R = z a (l- zfF(z), (10) 
then the radial equation satisfies the standard hypergeometric form |23| 

d 2 F dF 
z(l - z)-— + (c - (1 + a + &»-xr + a&F = 0, (11) 

where 

c = 2a + l, 
a6 = A + (a + /3)(a + (3 - 1), 
a + 6 = 2(a + /3)- 1 

and 

a 2 = -B, (12) 
/? = liVT - B. 

Without loss of generality, we put a = —iyB and (3=1 — \/\ — B. It is well known that the hypergeometric equation 
has three regular singular point at z = 0, z = 1 and z = oo, and it has two independent solutions in the neighborhoods 
of each point [2j| . The solutions of the radial equation reads as follows 



R{z) = Ciz a {\ - zfF(a, 6, c, z) + C 2 z- a {l - zfF(a, 6, c, z). 



(13) 
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In the gravitational case quasinormal modes of a scalar classical perturbation of black holes are defined as the 
solution the Klein-Gordon equation characterized by purely ingoing waves at the horizon $ ~ e -iu{t+r) ^ since at least 
classically outgoing flux is not allowed at the horizon. In addition, one has to impose boundary conditions on the 
solutions at the asymptotic region (infinity). It is crucial using the asymptotic geometry of the space time under 
study. In the case of asymptotically flat space time, the condition we need to imposes to the wave functions is a 
purely outgoing waves <& ~ e -^{t-r) a t the infinity j2J]. For non asymptotically flat space time (AdS space time for 
example), several boundary conditions have been discussed in the literature. In the two dimensional cases of BTZ 
black holes, the quasinormal modes for scalar perturbations were found in Refs. 21] 25] by imposing the vanishing 
Dirichlet condition at infinity: In Ref. |26j ] these modes were computed with the condition of a vanishing energy 
momentum flux density at asymptotia. In this paper we consider the vanishing flux condition. 

In the neighborhood of the horizon (z = 0) using the property F(a,b,d,0) = 1 the radial solution is given by 



R(z) = C x exp(-*£ ln(l - I± ) + C 2 exp(i^ ln(l - ^+ )), (14) 
la r la r 

while the boundary conditions of purely ingoing waves at the horizon demands C2 = 0. In order to implement boundary 
condition at infinity (z — 1), we use the linear transformation z — > 1 — z formula for the hypergeometric function and 
we obtain, 

r = c lZ a (i - 2) /3 r(c)r(c - a-b) a + b _ c + i _ + 

1 (c — a)l (c — 0) 

- zy-«-W T{c) ^\+ b - C) F(c - a, c - b, c - a b + 1, 1 - z). 
T(a)T(b) 

Using the condition that the flux is given by 



T 



2i 



{R*d lx R-Rd ll R*l 



(16) 



is vanished at infinity. Due to, B > 0, the asymptotic flux has a set of divergent terms, with the leading term of 
order (I - z) x - 2f) . Then according to Eq. (|I5(I each of these terms are proportional to 



T(c)T(a + b~c) 



T(a)T(b) 

and hence for a vanishing flux at z = I the following restrictions have to be applied 



(17) 



a = —n or b = —n, (18) 
where n — 0, 1, 2, .... These conditions lead directly to an exact determination of the quasinormal modes as follows 

i (n- l)(n + 3)a 

"-'2 ^TT ■ (19) 

Note that the quasinormal modes have the properties being purely imaginary and independent of the radius of the 
horizon of the sonic black holes. Besides this, it shows the instabilities of this kind of analog black hole under sonic 
perturbations. Due to, that pure imaginary frequency of the zero mode has the wrong sign (positive)that mean an 
exponentially growing mode. 

For the infinite limit damping (i.e. the n — * 00 limit) 



^00 = -i(n + ~)a 
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IV. CONCLUSIONS AND REMARKS 



In this paper we have computed the exact values of the quasinormal modes of Unruh's sonic black holes and 
according to Ref. || this QNM's should be behaved similarly to QNM's of near horizon Schwarzschild black holes. 
The QNM's are pure imaginary, this kind of QNM's was reported in Refs. H^|2^|2^. Therefore, for black hole that 
show this kind of QNM's is impossible to compute the area spectrum from the Kunstatter adiabatic invariant. On the 
other hand, the QNM's formula show that the zero mode exhibit an exponentially growing and therefore the Unruh's 
black hole in the limit that the background fluid smoothly exceeding the velocity of sound and, for the first order in 
the expansion becomes one instable geometry under sonic perturbations that excite the zero In opposite direction 
this result shows the large stabilities of the dumb hole for perturbations that excite the overtone with n > 1 and for 
all overtones in higher damping limit. We also want to note that the frequencies of the QNM's do not depend of 
radius of the horizon and on the angular momentum of the perturbation. They only have a proportional dependence 
the control parameter (a) which describes the velocity of the sonic flow, in according to Ref. [2(| where the decay 
rate for this geometry and the thermal emission was only proportional to the control parameter. 
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